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ABSTRACT 

Context. So far, multiple stellar systems harbor more than 130 extra solar planets. Dynamical simulations show that the outcome 
of planetary formation process can lead to various planetary architecture (i.e. location, size, mass and water content) when the star 
system is single or double. 

Aims. In the late phase of planetary formation, when embryo-sized objects dominate the inner region of the system, asteroids are also 
present and can provide additional material for objects inside the habitable zone (hereafter HZ). In this study, we make a comparison 
of several binary star systems and their efficiency to move icy asteroids from beyond the snow-line into orbits crossing the HZ. 
Methods. We modeled a belt of 10000 asteroids (remnants from the late phase of planetary formation process) beyond the snow-line. 
The planetesimals are placed randomly around the primary star and move under the gravitational influence of the two stars and a gas 
giant. As the planetesimals do not interact with each other, we divided the belt into 100 subrings which were separately integrated. In 
this statistical study, several double star configurations with a G-type star as primary are investigated. 

Results. Our results show that small bodies also participate in bearing a non-negligible amount of water to the HZ. The proximity 
of a companion moving on an eccentric orbit increases the flux of asteroids to the HZ, which could result into a more efficient water 
transport on a short timescale, causing a heavy bombardment. In contrast to asteroids moving under the gravitational perturbations of 
one G-type star and a gas giant, we show that the presence of a companion star can not only favor a faster depletion of our disk of 
planetesimals but can also bring 4-5 times more water into the whole HZ. 

Key words. Celestial mechanics - Methods: statistical - Minor planets, asteroids: general - binaries: general 


1. Introduction 

Nearly 130 extra solar planets in double and multiple star 
systems have been discovered to date. Roughly one quarter of 
these planets is orbiting close to or even crossing their system’s 
habitable zone (HZ), i.e. the region where an Ea rth-analog could 
retain liquid water on its surface llReinl 1201 4ll . While most of 
these planets are gas giants, the incredible ratio of one in four 
planets to be at least partly in the HZ seems to make binary star 
systems promising targets in the search of a second Earth, espe¬ 
cially for the next generation of photometry missions CH EOPS, 
TESS and PLATO-2.0. About 80 percent iIr^ IMl of the 
currently known plane ts in double sta r systems are in so-called 
S-type configurations llDvorald Il984l] . i.e. the two stars are so 
far apart so that the planet orbits only one stellar component 
without being destabilized. As most of the wide binary systems 
host not only one gas giant, their dynamical evolution is quite 
complex. The question whether habitable worlds can actually 
exist in such environments is, therefore, not a trivial one. 
Previous works on early stages of planetary formation show 
that planetesi mals accretion can be more difficu lt than in single 
star systems llThebault and Haghighinouil I20l4 and references 
therein]. This in turn can question the possibility of forming 
embryos in such systems. However, studies of late stages of 
planetary formation show that, should embryos manage to form 


despite these adverse conditions, then the dynamical influence 
of companion st ar s is no t prohibitive to form Earth-lik e planet s 
llRavmond et al.L l2004 iHaghighinour and Raymond boOTll . 
Eurthermore, it was shown that binary star systems in the 
vicinity of the Solar System are capable of s ustain i ng habitable 
worlds, once they are formed llEggl et all 120131 iJaime et^ 
l2014l] . As the amount of water on a planet’s su rface seems to 
be cr ucial to sustain a temperate environment IlKasting et al.L 
Il993l] . it is important to identify possible sources. Eor Earth, two 
mechanisms seem to be important: i) endogenous outgassing of 
primitive material and ii) exogenous impact by asteroids and 
comets sources have been established. Since neither can explain 
the amount and isotope composition of Earth’s oceans in itself, 
however, models that favor a combin ation of both sources seem 
to be more successful lllzidoro et ahll . The amount of primordial 
water that is collected during formation phases of planets in 
S-type orbits in binary st ar systems containing additional ga s 
giants has been studied bv IHaghighinour and Ravmondl ll2007ll . 
They have shown that the planets formed in a circumstellar 
HZ may have collected between 4 and 40 Earth oceans from 
planetary embryos, but main trend to appear was: the more 
eccentric is the orbit of the binary, the more eccentricity is also 
injected into the gas giant’s orbit. This in t u rn lea ds to fewer 
and dryer terrestrial planets. lOuintana et alJ ll2007ll emphasize 
that during the late stages of planetary formation, without the 
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Table 1. Physical properties of the secondary star in the binary 
system. 


Stellar-type 

M* [Mo] 

L* [Lo] 

T* [K] 

G 

1.0 

1.0 

5780 

K 

0.7 

0.38 

5200 

M 

0.4 

0.08 

3800 


presence of gas giant planets, binary stars with periastron >10 
an have a minimal effect on terrestrial planet formation within 
~2 an of the primary, whereas binary stars with periastron < 
5 an restrict terre strial planet formation t o with in ~1 an of 
the primary star. lOuintana and Lissaued ll2006ll studied the 
late stages of planetary formation in P-type orbits in binary 
star systems (with maximum separation equal to 0.4 au and 
maximum secondary’s eccentricity equal to 0.8) including 
Jupiter and Saturn-like planets. They conclude that the higher 
the secondary’s apoapsis, the smaller the number of planets 
formed and the lower their mass. The anti-correlation between 
a system’s eccentricity and the number of planets has also been 
fo und in a preliminary interpr etation of observation statistics 
by iLimbach and Turnej 1201411 . These results indicate that the 
likelihood of finding habitable planets in such an environment 
could be small. 

Stochastic simulations proved that planet can also be formed 
almost dry in the circumprimary H Z of binary star systems 
iHaghighipour and Raymond l2007ll . However, as emphasized 
by these authors, water delivery in the inner solar system is not 
only due to radial mixing of pla netary embryos. Smalle r objects 
can also contribute as shown in iRavmond et al.l ll2007ll . Indeed, 
they highlighted that water delivery from smaller planetesimals 
is statistically robust and should supply terrestrial planets with 
a significant water source of perhaps three to 10 oceans. In our 
work, we aim to answer how much water can be transported into 
the HZ via small bodies, remnants from the late phase of the 
planetary formation, providing thus other water sources to em¬ 
bryos. We statistically study the dynamics of an asteroid belt in 
such systems and we treat this problem in a self-consistent man¬ 
ner as all the gravitational interactions in the system as well as 
water loss of the planetesimals due to outgassing are accounted 
for. Our main goal is to examine the influence of the secondary 
star on the flux of small bodies from icy regions to the HZ and 
the efficiency of the water transport within a short timescale of 
lOMyrs. 


2. Initial conditions and dynamical model 

Our study is focused on a primary G-type star with mass equal 
to one solar mass and we investigate the dynamical effect of 
a secondary of either G, K or M-type with physical properties 
expressed in Table[T] The studied binary star systems encompass 
relatively tight configurations, i.e. semi-major axes in a range 
of ab e [25 : 100] au. This parameter has been changed in steps 
of 25 au in our simulations. The secondary is on an elliptical 
coplanar orbit with eccentricities eb e [0.1 : 0.5] increased in 
steps of 0.2. In total 12 configurations have been investigated for 
any given stellar type of the companion. In addition, we consider 
a gas giant located at aco = 5.2 au on a circular coplanar orbit, 
with a mass equal to Jupiter’s mass. 


Table 2. Mass loss in the simulated collisions. Mass loss is de¬ 
fined as the quantity not in the surviving bodies (one for erosion, 
two for hit-and-run and merging) after one impact. 


Collision scenario Mass 


V 

[^esc] 

a 

n 

Type 

[wt-%] 

1.00 

20 

merging 

0.6 

1.27 

33 

merging 

0.9 

1.37 

42 

hit & run 

0.4 

1.99 

44 

hit & run 

3.5 

2.77 

22 

erosion 

75.2 


A disk of planetesimals is modeled as a ring of 10000 asteroids 
with masses similar to main belt objects in the Solar System and 
each asteroid was as signed an initial water mass fraction (here - 
after wmf) of 10% llAbe et al.L l2000i iMorbidelli et al.L l2000ll . 
During the late phase of planetary formation, if the inner region 
of the system is mainly dominated by large embryos (follow¬ 
ing a speci fic mass distribution) with masses between Moon to 
Mars size llRavmond et al.L l2004 iHaghighinour and Raymond 
|2007|] . debris resulting from collisions of such embryos are 
also present. To determine the lower and upper limits for 
the asteroids’ masses, we performed independent preliminary 
simul ations wi t h a 3D smooth par t icle hy drodynamics (SPH) 
code llSchaferl 120051 iMaindl et al.L |20T3D. First-order consis- 
tency is ach i eved b y a tensorial correction as discussed in 
ISchafer et al.l Il2007ll . It includes self-gravity and models mate¬ 
rial strength using the full elasto-plastic continuum mechanics 
and the Grady-Kipp fragmentat i on model for fracture and brittle 
failure IlGradv and Kippl Il980l iBenz and Asphau^ Il994|] . The 
scenarios involve collisions of rocky basaltic objects with one 
lunar mass at different encounter velocities and angles a. The 
latter are defined in a way so that a = 0° corresponds to a head- 
on collision. The scenarios start with the bodies five diameters 
apart to let the SPH particle distribution settle. This prelimi¬ 
nary simulation time-span was 2000 min. In an earlier study, 
most collisions of Moon-sized bodies in the Solar System’s 
HZ were found to happen at velocit i es v $ 2vesc at arbitrary 
collision angles IlMaindl and Dvor^ |2014|] . We chose initial 
conditions in this range which are in the merging/partial ac- 
creti on and hit-and-run domain of the collision outcome map 
[cf. Agnor and A sphaugl |2004 iLeinhardt and Steward 1201 2l 
iMaindl et al.Ll2014ll . Expecting a somewhat higher spread in v in 
binary systems we also included a scenario in the erosion/mutual 
destruction domain. Table lists the collision scenario parame¬ 
ters and gives the resulting mass loss of the survivors (one body 
for merging, two bodies in the erosion and hit-and-run scenarios, 
respectively) after one impact. 

The fragment sizes beyond the survivors drop significantly in 
the hit-and-run and merging scenarios (Fig.[T^, b): In the mutual 
destruction case, all of the ten largest fragments possess masses 
^ 1 % of the total system mass, which is ~ Ceres’ masfl and 
hundreds are above th e “sign ificant” fragment threshold in the 
sense of IMaindl et akl Il20l4 . The smallest fragment consists 
of one SPH particle (0.001% of the total mass for 100k SPH 
particles) which corresponds to ~ 0.1% of Ceres’ mass. As 
increasing the number of SPH particles will result in even 
smaller fragments, this mass is an upper limit for the smallest 
fragment. As this latter will contain ~ 0.006% Earth-oceans 
unitfl we chose to neglect the water contribution of smaller 


* Ceres’ mass is equal to 4.73x10 ‘“Mg 
^ 1 ocean = 1.5 x 10’' kg of HjO 
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Fig. 1. Collision scenarios resulting in a merge (a), a hit-and-run encounter (b) and erosion/mutual destruction (c). The diagrams 
show the masses of the ten largest fragments at the end of the simulation as percentages of the total system mass Mtot = 2 lunar 
masses. The dashed horizontal line indicates the limit for significant fragments (see text). 



particles. Our minimum and maximum mass are thus defined 
according to the fragments’ mass after one impact. Therefore, 
members of our ring will have masses randomly and equally 
distributed between 0.1% to one Ceres’ mass. The total mass of 
the ring of 10000 asteroids amounts to Mr =0.5 M@. Thus the 
quantity of water, in terms of Earth-ocean units, available in a 
ring will be 200. 


In each system, we defined the size of the disk of planetesimals 
by the following borders. As we focus on the transport of icy 
bodi es, the inner border of the disk is set to the snow-li ne posi¬ 
tion llLecar et al.L 120061 iMartin and Livioi 1201 2i 1201 3ll . border 
between icy and rocky planetesimals. Its position changes as 
the star’s luminosity evolves during its birth phases. As our host 
star is a G2V type, we considered the inner border of the disk 
of planetesimals according to observations in the Solar System, 
suggesting the snow-line to be at ~ 2.7 au. The outer border of 
the disk of planetesimals is influenced by perturbations induced 
by the secondary star. Since we consider initially circular 
motion for the planetesimals, the stability border depends 


mainly on three paramete r s: the mass ratio of the sys t em, ai, 
and eh llRabl and Dvorald I1988L iHolman and Wiegeri. 1 19991 
iPilat-Lohinger and Dvorakl 2002ll . These authors showed 
that it is possible to link these parameters to derive a critical 
semi-major axis a^ as the maximum initial se mi-major axis for a 
particl e to su rvive in the system. In contra st to lRabl and Dvor^ 
ill988ll and iHolman and Wieg^ lll999ll who classified un- 
stable orbits via ejections of te st planets from the system, 
iPilat-Lohinger and Dvorakl ll2002ll calculated the Fast Lyapunov 
Indicator (ELI) for each orbit to distinguish between stable 
and chaotic motion . This well known cha os detection method 
was introduced by iFroeschle et al.l lll997ll . In case of circular 
motion of the pl anets and the planetesimals , it is possible to 
use the study by IHolman and Wieg^ Ill999l] where ac and its 
uncertainty Aaj. allows a good determination of the outer border 
(maximum semi-major axis allowed) for asteroids in the ring 
as a^ - Aac, which is in good agree ment with the stability limits 
derive d from ELI computations llPilat-Lohinger and Dvorakl 
I 2 OO 2 II . Asteroids of our belt are randomly positioned between 
the inner (the snow-line) and and outer borders (the stability 
limit a^ - Aa^). Figure |2] shows a comparison of initial asteroid’s 
mass distribution for a secondary M star at 100 au and Cb = 0.1 
and 0.5. The position of the gas giant planet is indicated by 
the vertical line. Because of the binary’s tight periapsis in the 
case of eb = 0.5, the asteroid belt is less extended and exhibits a 
higher mass density. In order to prevent immediate dynamical 
instability, all bodies have been carefully placed so that their 



a [au] 


Fig. 2. Example of mass distributions in the circumprimary ring 
as a function of the asteroids’ semi-major axis under the gravi¬ 
tational influence of a secondary M star at ab = 100 au and Cb = 
0.1 (■) and Cb = 0.5 (•). The vertical line refers to the gas giant’s 
position. All particles are initially spaced by several Hill’s radii 
to prevent immediate dynamical instability. 


initial mutual separation were at least several Hill’s radii. 

In our main simulations, as we assume that the giant planet is 
already formed so that the remaining gas has been evaporated or 
coagulated into the asteroid belt, we do not take any effects re¬ 
lated to gas drag into account (therefore no gas driven migration, 
no eccentricity dampening). The initial eccentricities and incli¬ 
nations are randomly chosen below 0.01 and 1° respectively. To 
avoid strong initial interactions with the gas giant, we assume 
that the giant planet has cleared a path in the disk around its or¬ 
bit. The width of this gap is +3Rh,gg where Rh.gg is the giant 
planet’s Hill Radiu]^- Requiring orbital stability of the gas giant 
at 5.2 au reduces the number of possible binary configurations. 
Therefore, we excluded the case where ac - Aac < aoG + 3Rh,gg- 
This results in a total number of 23 binary systems configura¬ 
tions that were studied in this work and summarized in Table |3] 

We limited our study to 10 Myrs integration time. This is of 
course much shorter than the timescale of terrestrial plane- 

^ However, we did not exclude possible location of asteroids inside 
mean motion resonances (MMRs). Indeed, we assume that the presence 
of the gas might have kept asteroids on stable orbit inside MMRs. Then, 
the MMR perturbations became stronger when the gas vanished. 
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Table 3. Binary configurations studied in this article. Only the 
secondary stellar-type is reported. 


e/, 

25 

50 

75 

100 

0.1 

K-M 

G-K-M 

G-K-M 

G-K-M 

0.3 


G-K-M 

G-K-M 

G-K-M 

0.5 



M 

K-M 


tary formation (~ 100 Myrs) but as we make a comparison 
of many different binary star systems, we had to restrict this 
study to a shorter time. However, we will select from this 
study the most interesting systems for which a statistic over 
100 Myrs will be made. W e numerically integrated our sys¬ 
tems using the nine package llEggl and Dvorakl 1201 011 and only 
gravitational perturbations were taken into account. The nu¬ 
merical integrat or used for the computa ti ons is based on Lie- 
series (see e.g. iHansl meier and Dvorakl lll984ll and more re- 
centlv lBancelin et al.l II 2 OI 2 I] '). For a given configuration in Table 
[3] as our planetesimals do not interact with each other, the disk 
was divided into 100 subrings and separately integrated. 


3. The Habitable Zone borders 


As we study the flux of water-rich asteroids into the HZ, we 
need to know the position of this area. The definition, modeling 
an d compu t ation o f the ’’clas s ical H Z ” (hereafter CHZ) is g iven 


in 


___ Kasting ill988 1. K astingl lll99lll . iKasting et al.l 11199311 and 
iKoDParapu et al.l I 2013ll . All these studies are based on a ID 
cloud-free climate model, where the inner edge of the HZ is 
computed by increasing the surface temperature and the outer 
edge, by increasing the CO 2 partial pressure (maintaining a 
constant surface temperature at 273 K). The corresponding 
stellar flux, needed to maintain the surface temperature, is then 
derived. Thus, these authors were able to express the CHZ 
borders as a simple fit func tion containing the stell ar luminosity 
and temperature. Recently. iKopparapu et al.l ll2014ll investigated 
the dependence of the HZ borders on the planetary masses and 
derived new coefficients for the computations of the effective 
solar flux. For a G star, the inner edge of the HZ corresponding 
to the runaway greenhouse limit is 0.950 au and the outer edge 
of the HZ corresponding to the maximum greenhouse is 1.676 
au. The updated i nner edge value is clo ser to the Sun than 
the one found by iKopparapu e^alj j 2Q13ll b ecause they used 
inputs from a 3D model by Ikeconte et akl Il2013l] . However, 
the exact value for the inner CHZ border depends on many 
assu mptions and is still u nder discussion in current literature, 
fe.g. lWolf and Toonl Iml . 


When including a binary companion in the system, the ad¬ 
ditional gravitational interaction an d radiation can shrink 
the CHZ borders, as shown in lEggl et all Il2012ll and 
iKaltenegger and Haghighipouri ll2013ll . As a matter of fact, 
the smaller the periapsis of the binary, the more important the 
insolation of the primary, as the periapsis distance of the planet 
with respect to its host st ar can change sig nificantly. In order 
to account for this effect, lEggl et aP ll2012ll introduced the so- 
called Permanently Habitable Zone (hereafter PHZ). The PHZ 
contains information on the planet’s perturbed orbit. It serves as 
a mean to distinguish areas where a planet always receives the 
correct amount of insolation to remain habitable from regions 
where insolation conditions for habitability are only fulfilled 
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Fig. 3. PHZ/CHZ ratio as a function of the binary separation 
and the secondary’s eccentricity. The red color indicates that the 
planet will always lie outside the PHZ borders. Only a graph 
for a secondary G star is represented as it is quite similar for a 
secondary K and M-types. 


in an avera sense. Figure |3] shows the ratio PHZ/CHZ as a 
function of the binary’s orbital elements (ab,eb) for a G2V-G2V 
binary star system. The blue and black colors cover a large 
region where no or only minor differences between PHZ and 
CHZ were found. This means that the additional radiation from 
the second star is not enough to drastically cause a change of 
insolation at the surface of a planet. This is due to the fact that 
the secondary is too far or its periapsis too far away. For higher 
eccentricities, where the secondary’s orbit approaches closer 
to the host star, the difference is increased (green color). The 
red color refers to regions where the PHZ vanishes either due 
to excessive insolation or due to orbital instability. It is clear 
that the truncation of the PHZ increases for large values of the 
periapsis of the binary companion. This is due to the fact that 
the stellar gravitational perturbations acting on the planet will 
cause a significant change of the planet’s periapsis, which in 
turn influences the insolation at the planet’s surface. A similar 
behavior is observed for K and M class secondaries as their 
mass does not significantly influence the PHZ/CHZ ratio. For 
our binary configurations, the difference between CHZ and PHZ 
is only of the order of 10% as shown in TableH) As these values 
do not vary strongly with the secondary’s mass, we did not 
report all 23 configurations. 


4. Asteroid flux and water transport to the HZ 

4.1. Statistics on the disk dynamics 

During the simulation, each particle is tracked until the end of 
the integration time in order to assess the following numbers; 

• asteroids crossing the HZ. They will be referred to as 
Habitable Zone crosser^ (hereafter HZc). As we assume a 
two dimensional HZ, an asteroid will be considered as a HZc 
if the intersection point between its orbit and the HZ plane 
lies within the HZ borders. 

^ The planet’s orbit is eccentric enough could leave the HZ from time 
to time. 

^ A HZc can cross several times the HZ before leaving the system or 
colliding with the stars or the planets 
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Table 4. PHZ borders given in an as a function of the binary orbital characteristics (at, Cb) used in our simulations. Compared to the 
CHZ borders (inner border is 0.950 au and outer border is 1.676 au), the PHZ intervals are smaller due to the additional radiation 
from the second star and eccentricity injected in the planet’s orbit. 


[au] 

25 

50 

75 

100 

0.1 

[0.959:1.654] 

[0.955:1.664] 

[0.954:1.668] 

[0.953:1.670] 

0.3 


[0.964:1.639] 

[0.959:1.651] 

[0.957:1.657] 

0.5 



[0.968:1.627] 

[0.964:1.638] 


• asteroids leaving the system when their semi-major axis> 
500 au 

• asteroids colliding with the gas giant or the stars 

• asteroids still alive in the belt after 10 Myr 

Figure |4] shows the resulting statistics on the asteroids’ dynam¬ 
ics. It represents four panels corresponding to the secondary’s 
semi-major axis investigated. Each histogram shows the dynam¬ 
ical outcome of our asteroids expressed in terms of probability, 
as a function of the secondary’s eccentricity. Below 100%, the 
percentage of asteroids that are still present in the belt (’’alive”), 
that were ejected or collided with the stars or the gas giant is 
shown. The black area of each histogram above 100% indicates 
the probability for asteroids to enter the HZ. Such asteroids, 
crossing the HZ, are called HZc. The perturbations due to both 
the gas giant and the binary companion can cause an increase in 
eccentricity of the asteroids within the planetesimal disk which 
of course depends on the binary configuration. A comparison 
of the different histograms indicates that the key parameter 
is the periapse of the binary system which is defined by the 
eccentricity of the binary. Of course not only the asteroids 
are perturbed by the secondary, the gas giant at 5.2 au is also 
perturbed and its initially circular motion will change to an 
elliptic one. This behavior is highlighted by the increasing 
value of the probability to become a HZc if the secondary’s 
periapsis distance decreases and if its mass increases. Indeed, 
for a given value of a* (for instance 50 au), one can see that 
this probability is at least doubled when e* increases. As a 
consequence, the asteroid belt will be depopulated because of 
dynamically induced ejections as well as collisions with the 
giant planet and the stars. Because the rate of colliding and 
ejected asteroids is increased, a ring will be depopulated faster 
when ej becomes larger. Therefore, the statistics in Fig.|4]shows 
that the probability for a member of the asteroid ring to stay 
in the system after 10 Myrs will decrease with the periapsis 
distance and the secondary star’s mass. Finally, these results 
can also answer the question of the presence of an asteroid 
belt in such systems. Indeed, if we assume that the gas could 
protect small bodies from secular resonances or mean motion 
resonances (MMRs), it is highly unlikely that they survived 
in binaries with small periapsis separation, after the gas has 
dissipated. 


4.2. Timescale statistics 

Depending on the periapsis distance of the secondary, the disk of 
planetesimals can be perturbed more or less rapidly. As a matter 
of fact, asteroids will suffer from the gravitational perturbations 
of the secondary star and the gas giant, and their eccentricity 
may increase quickly. Figure |5] shows statistical results of the 
average time needed by an asteroid to become a HZc, i.e. the 


eh = 0.1 

eh = 0.3 


75 

lb [au] 


Fig. 5. Median time for an asteroid to become a HZc for the case 
of a secondary G-type with Cb = 0.1 (■) and Cb = 0.3 (•). This 
corresponds to the time when an asteroid crosses the HZ for the 
first time. The statistics is made over the 10000 asteroids and the 
Icr value is represented by error bars. 


time it takes to reach the HZ. This corresponds to the time spans 
until the first asteroid enters the HZ. The median value and its 
absolute deviation (error bars) are presented for a set of 10000 
asteroids for the case of a secondary G-type with eb = 0.1 (■) 
and eb = 0.3 (•). This confirms a strong correlation between the 
periapsis distance and the time of first crossing. Figure|5]clearly 
shows that the average time varies from a few centuries to tens of 
thousands of years. The closer and more massive the secondary 
star, the sooner asteroids can reach the HZ. In contrast, the cross¬ 
ing timescale will become larger as the number of HZc increases 
as shown in Fig.|6l This parameter corresponds to the bombard¬ 
ment timescale, within 10 Myr of integration time, and is derived 
regarding the last crossing inside the HZ but regardless the as¬ 
teroid’s water content. Indeed, once an asteroid becomes a HZc, 
it can cross the HZ several times as long as its orbit is stable 
and until it is ejected out of the system or collides with the giant 
planet or the stars. However, water-rich asteroids could be dry 
before this corresponding timescale, revealing that water trans¬ 
port could occur on a very short time (compared to planetary 
formation timescale). Besides, one can notice from Fig. |6] that 
most of the systems, with large a* and low e*, having low cross¬ 
ing timescales are also those where most asteroids remained in 
the after 10 Myrs as shown in Fig. 01 This suggests that asteroid 
flux to the HZ could occur in several steps as some asteroids need 
more time to drastically increase their eccentricity in order to be 
moved to lower orbits and reach the HZ. This reveals also that 
water sources can still be available in the ring, provided that as¬ 
teroids still have icy water on their surface. This study would of 
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Fig. 4. Statistics on the disk of planetesimals dynamics. Each histogram shows the evolution (expressed in probability) of the 
asteroids in the ring within 10 Myr of integration. They can still be present in the system (’’alive”), become a HZc, collide with the 
stars or the gas giant, or be ejected out of the system. The closer the secondary and the higher its mass, the higher the probability to 
empty the asteroids ring. 


course require longer integration time > 10 Myrs. Nevertheless, 
we have a clear guess on the efficiency of binary stars systems 
to transport asteroids from beyond the snow-line to the HZ on a 
short timescale. 

4.3. Water transport statistics 

Each HZc entering the HZ will bring a certain amount of water. 
However, regarding the integration timescale, the water content 
of asteroids may vary with time. Indeed, increased eccentricities 
can cause asteroids to approach close to the stars. This in turn 
would lead to a mass variation mainly due to a loss of their water 
content. 


4.3.1. The water mass loss process 


To quantitatively assess the water content of asteroids, we fol¬ 
lowed their water mass fraction evolution throughout the simu¬ 
lations, including mass loss process, until they enter the HZ. The 
main mechanisms that can induce a relevant mass loss for active 
(comet-li ke) or i nactiv e asteroids are ice sublimation and impact 
ejection llJewittl l2012ll . The latter process occurs when smaller 
asteroids impact larger onefl These impact s can be highly ero- 
sive due to characteristic speeds ~ 5 km.s ‘ llBottke et al.Lll994|] 
and the ejected amount mass can be significant. However, typ¬ 
ical impact events in t he main-belt have an impa ct probability 
~ 3.0 X 10~^ km~^.yr~‘ llEarinella and Da^ll992|] . Yet, accord¬ 
ing to lBottke et al.l Il2005ll . the timescale for an impact event to 


happen in our sample of asteroids belt (with radius from tens to 
hundreds of kilometers) is much longer that our integration time. 
Thus, we neglected this process in our study. The only mass loss 


® Such impacts do not necessarily lead to a complete destruction or 
break-up of asteroids 


process we consider is due to ice sublimation when the asteroid 
comes relatively close to the star. This process can be stepped up 
in double star systems, especially when the companion is on an 
eccentric orbit. Even if the secondary is not particularly close, 
asteroids’ eccentricity can be pumped and they can receive a 
large amount of insolation. Therefore, the water mass loss rate ih 
was computed accounting for the radiation of both stars. The es¬ 
timation of lit can be derived by solving a balance energy equa¬ 
tion between the net incoming stellar flux, the ice sublimation 
and the thermal re-radiatiorQ, as expressed in Eq. ([1]) 

^E*i(l-A.) 

> - - cos 6»i = eo-T'* -H L(T) m(T) (1) 

4-J R2[au] 


where 

• E* is the stellar constant and is computed as E*=EoL* (Eo 
~ 1360 W.m ^ is the solar constant). 

• A, the bound albedo, product of the geometric albedo and the 
phase integral, defines the fraction of the total incident stellar 
radiation reflected by an object back to space. Asteroids can 
have 0.5 but most of them h ave relatively low albedo 
BShestopalov and Golubevall201 111 . Ice material is known to 
be a good radiation reflecto r. In order to maximize m, we 
will follow the approach of iJewitj Il2012ll considering our 
objects with dirty ice material with low albedo because clean 
ice sublimates too slowly at main-belt distances. Therefore, 
we used an averaged albedo A, -A 2 -A- 0.05. 

• R is the distance to the star (primary or secondary) expressed 
in au 

• 0 is the angle between the incident ligh t and the norm al to 
the the asteroid’s surface. According to iJewitj Il2012ll . m is 
weaker for an isothermal surface than at the sub-solar point 

^ assuming the asteroid to have black body properties 
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Fig. 6. Crossing timescale expressed in Myr. It corresponds to the last crossing inside the HZ. Up to 6 Myr for the most perturbed 
systems but also with the lowest number of asteroids still present in the system within 10 Myr of integration time. For the other 
systems, the crossing timescale is much lower because their belt still contain a huge number of ’’alive” asteroids. Indeed, the 
dynamics of the system needs more time to empty the asteroid belt population. 


of a non-rotating body. Again, to maximize rii, we consider 
the sub-solar case with 0 = 0°. 

• e is the emissivity of the surface, e ~ 0.9 

• cr = 5.67 X 10'* W.m'^.K'^ is the Boltzmann constant 

• T is the equilibrium temperature at the surface expressed in 
K 

• L is the latent heat of sublimation in J.kg'* 

• ih is the surface mass loss rate in kg.m'^.s'* 

m can be obtained by computing the temperature T solving Eq. 
O. To this purpose, we expres sed all the variable parameter s as 
a function of T. According to iDelsemme and Millet! Ill97lll . m 
can be written as: 

where // = 18 g.mol'* is the water molar mass and Ps is the sat¬ 
uration pressure and defined by the empirical formula: 

2484.986 

logPs = 4.07023- - --H 3.56654 log(T) -0.00320981 T 

(3) 

Finally, the latent heat is expressed as: 

L (T) = 2834.1 - 0.29 (T- 273.15) - 0.004 (T- 273.15) J.g'* 

(4) 


4.3.2. The water mass fraction of incoming HZc 

The value of m is constantly updated during the simulations and 
the accumulated surface mass loss dm reads: 

dm-'^mAt kg.m'^ (5) 

where At represents the time elapsed since the beginning of the 
integration. We then compute the sublimating area 47rr^ with r 


the radius of the HZc. If we consider the following density for 
water ice shell and a basalt core, respectively p, - 900 kg.m'* 
Pc = 3000 kg.m'*, then for a given wmf, the mean density of an 
asteroid is 

p = (1 - wmf) Pc -H wmf X Pi 

Therefore 



with m the mass of the asteroid. Thus we can derive the total 
water mass loss in kg Am = dmxAnr. 

Whenever an asteroid becomes a HZc i.e. when it first enters the 
HZ, we suppose that its current water content will be delivered. 
This is the best case scenario as the total water content of the 
HZc corresponds to the maximum amount of water available in 
the HZ. In reality, however, only a small fraction of water would 
be accreted onto a planet (or planetary system) inside the HZ. 
Indeed, both the real position of the planet on its orbit and the 
impact velocity of these HZc are ke y parameters in collis ional 
material transport as pointed out by iThebault et alJ Il2006l] and 
iLeinhardt and Stewart! Il2012ll . A fully self consistent modeling 
of the water delivery at impact lies beyond the scope of the 
current work, however. Hence, we aim to study the amount 
of water transported into the HZ rather than the amount of 
water accreted by a planet therein. Figure |7] shows the total 
amount of water that was transported into the HZ (expressed 
in terrestrial oceans units) as a function of the binary system 
characteristics. Each histogram represents the total amount of 
oceans ending in four partitions of the HZ. They correspond 
to equally spaced rings using the values obtained in Table |4] 
Therefore we define the inner HZ border (closest ring to the 
PHZ inner border value). Central 1 and 2 (intermediate borders) 
and Outer HZ (closest ring to the PHZ outer border value). We 
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also illustrated on this figure, the additional quantity of water 
transported between the inner edge A [CHZ{Ay,PHZ(A)] and 
the outer edge B [PHZ{B)\CHZ{B)] of the HZ (gray color and 
refereed as ’’additional”)- 

The main results suggested by this figure are; 


a) one can see that a maximum of 70 oceans can be transported 
into the CHZ. Note that 70 oceans represents 35% of the 
total number of oceans contained in the 10000 asteroids 
populatioifl Variances are large, however. One of the 
reasons for the large spread in the amount of transported 
water is the short simulation time. This introduced a bias 
towards systems where the water transport is fast. In other 
words, the seemingly small amount of transported water 
in some of the presented systems does not imply that their 
HZs have to remain dry. It merely highlights the fact that 
the delivery process would require a longer time as many 
asteroids are still present in the system after 10 Myrs of 
integration time (see Figure Ell. However, these results show 
how fast and efficient are some systems to transport water. 

b) the difference between the water transported in CHZ and 
PHZ (gray color) is not that significant compared to the total 
number of incoming oceans. Indeed, the difference does 
not exceed ~3 oceans. One should note that we disregarded 
large orbital variations of any embryos located in the HZ 
caused by the perturbed motion of the gas giant because of 
the presence of the secondary star. This in turn can perturb 
the motion on any plane t s in th e H20. However, according 
to IWilliams and Poliak Il2002|] . an average insolation - 
covering almost the entire CHZ - is sufficient to retain liquid 
water on an Earth-twin surface and secure the habitability of 
planets. 


c) any binary star systems is efficient enough to produce a flux 
of asteroid within the whole HZ and thus to make water 
sources available for embryos lying there. However, we can 
see that statistically, most of the water ends in the outer 
HZ. This is due to the fact that its surface area is much 
more wider than the other rings and the probability for an 
asteroid to fell in the outer HZ is higher. One can notice 
that the quantity of water brought into the outer HZ ver¬ 
sus the other cells is quite balance. Indeed, regarding the 
definition of a HZc, we expect asteroids on inclined or¬ 
bits not to necessarily cross the outer HZ. Besides, dynam¬ 
ical studies suggest that secular perturbations can lie inside 
the HZ or beyond the s now-line, depending on the binary 
systems’ characteristics iBancelin et al.L 1201 51 iBaszo et’aU 


l2015Llraat-Lohinger et all 201511 . In addition, depending on 

the secondary’s periapsis, the seco ndary will shorten th e 
lifetime of particles inside MMRs llBancelin et all 1201 5ll . 
Therefore, for asteroids initially orbiting inside MMRs 
and/or secular perturbations, the dynamical evolution can be 
violent and their orbit crossing the HZ will not necessarily 
be from the outer HZ to the inner HZ. 


° 200 oceans contained in this population 

® when including a gas giant in a binary star systems, its additional 
perturbation will increase the eccentricity on any planets located in the 


HZ 


4.4. Influence of the initial water content 

Let us now compare the water transport efficiency when 
the initial wmf is not equally distributed throughout the 
asteroid belt. The equal water distribution model will be 
called WMFa- Observations in the main-belt llAbe et all 120001 
iDeMeo and Carrvl 1201 4|] suggest that a gradient exists in the 
chemical composition of the asteroids. Besides, we expect dis¬ 
tant asteroids, up to the Kuiper-belt distance, to have a higher 
wmf than asteroids in the main-belt. Thus, we modeled - we 
called this model WMFb - the wmf distribution as a linear 
function of the distance to the central star, fulfilling the fol¬ 
lowing boundary conditions: the upper limit is fixed at 20% 
of water for asteroids at 30 at r^ To find the l ower limit 
for asteroids at 2.7 au, we f ollo''^Ravmond et alJ Il2004|] and 
iHaghighipour and Raymond Il2007l] assuming a wmf ~ 5%. 
Under these assumptions, a belt in tight binary systems will have 
a lower water mass ratio than a belt under the influence of a sec¬ 
ondary on a low eccentric orbit. Indeed, the ratio between model 
WMFb and WMFa gives lower limits of 0.62, 0.58 and 0.63 re¬ 
spectively for G, K and M secondary stars. When the secondary’s 
periapsis increases, we get upper limits of 1.07, 1.18 and 1.27 
respectively. However, the quantity of oceans transported to the 
HZ using WMFb does not exceed 2/3 of the amount of trans¬ 
ported water when using WMFa, even if the belt initially con¬ 
tains more water. This shows that: 

a) our results are robust 

b) asteroids closer to the snow-line are more likely to become 
HZc. As a matter of fact, according to model WMFb aster¬ 
oids will contain 10 % of water if they are located below 11 
aiO If asteroids becoming HZc were initially beyond this 
distance, we would have had an equivalent or higher number 
of transported oceans. 


5. Comparison with a singie G star system 

In this section, we will compare the water transport efficiency 
between binary and single star systems. To this purpose, we con¬ 
sidered the same initial conditions for the gas giant and the as¬ 
teroid belt distribution, except that we only have a single G star 
in our dynamical system. As the comparisons are made for the 
same asteroids’ initial conditions, we have to consider, for the 
single star system, the same size for the disk, given by the bina¬ 
ries’ characteristics. This is why, in Fig. |9] results for the refer¬ 
ence single star case are different depending on the values of the 
binary’s eccentricity and separation. 

For computational reasons, our comparison is limited to a 
G2V-G2V-type binary. Besides, regarding Fig. |6] there is no 
more crossings in the HZ after 5-6 Myr. Therefore, the number 
of oceans brought to the HZ will not vary after this time. Thus, 
each system was compared at equivalent integration times (< 5 
Myr). Our results show that an asteroid that was initially in the 
ring will need 2-20 times longer to reach the HZ in a single 
star system. That is, because the asteroid belt is not perturbed 
strongly enough by the gas giant to produce a large asteroids 
flux towards the HZ as shown in Fig. [8j the presence of the sec¬ 
ondary star can put asteroids with initial semi-major axis ao, on 
high eccentric orbits, i.e. very low periapsis distance qnzc- This 

As shown in Figure [T] our objects’ semi-major axis does not go 
beyond 30 au 

provided that the critical semi-major axis criteria allows asteroids 
to be located beyond this distance 
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Fig. 7. Total number of oceans transported by all the HZc when they first cross the HZ. The color code is related to four equally 
spaced rings of the HZ: inner, central 1 and 2, outer HZ (see text). We also indicated the equivalent number of additional ocean 
crossing the truncated CHZ interval when the dynamics of the secondary star is not taken into account for the computation of the 
HZ borders. 
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Fig. 8. Initial semi-major axis ao and final periapsis distance qnzc 
of the incoming HZc for a single (•) and binary (■) star system. 


will result in a faster depletion of the belt and a shorter bombard¬ 
ment duration, compared to a single star system. 

Consequently, the probability for an asteroid to cross the HZ 
will be ~ 4 - 6 times higher in a binary star system. Note that 
the number of scattered incoming asteroids mainly increases be¬ 
cause the secondary star perturbs the giant planet’s orbit. As 
the flux of asteroids is less efficient in a single star system, the 
amount of water brought to the HZ is smaller. In fact, the wa¬ 
ter transport is 4 - 5 times less efficient without a second star. 
Finally, Fig.^compares the amount of delivered water in various 
systems. Histograms marked with the letter B refers to binary 
star systems, those with letter S to single star systems. The color 
code indicates the amount of water that ended up in four equally 
spaced sub-rings (Inner, Central 1, Central 2, Outer) of the corre¬ 



<!b =b 

Fig. 9. Comparison of the water transport in single star (S) and 
binary star (B) systems. We restricted this study to binaries with 
two G-type stars. The various panels correspond to different 
semi-major axis of the secondary: at = 50 - 75 - 100 au. The 
color code refers to different sub-rings of the HZ (see text). Note 
that, the two results for each panel account for different initial 
asteroid belt distributions (depending on Ob), taken as the same 
for (S) and (B). 


sponding PHZ (see Table H]). For a single star system, each ring 
of the CHZ is computed using the inner edge value 0.950 au and 
outer edge value 1.676 au. It is not surprising that the outer HZ 
is the most crossed ring. Indeed, its area is much larger than the 
other rings. This figure also highlights the fact that in such sin¬ 
gle star systems, hosting only one giant planet, basically all the 
water is transported in the outer HZ because the perturbation is 
not strong enough to drastically increase the eccentricity of any 
asteroid in the belt. In any case, these results show the efficiency 
of a binary star to transport water in the entire HZ over a shorter 
timescale compared to a single star system. 
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6. Summary and conclusion 

In this work, we investigated the influence of a secondary star 
on the flux of asteroids to the Habitable Zone (HZ) over 10 
Myrs of integration time. We estimated the quantity of water 
brought by asteroids located beyond the snow-line into the HZ 
of various double stars configuration (separation, eccentricity 
and mass of the secondary). An overlap of perturbations from 
the secondary and the giant planet in the primordial asteroid 
belt causes rapid and violent changes in the asteroids’ orbits. 
This leads to asteroids crossing the HZ soon after the gas in the 
system has dissipated and the gravitational dynamics become 
dominant. Our results point out that binary systems are more 
efflcient for water transport into the HZ, compared to a single 
star system. Not only asteroids flux is 4 - 6 times higher when 
a secondary star is present, but also the amount of transported 
oceans to the HZ can be 4 - 5 times more important, providing 
other water sources to embryos, in the whole HZ, in the late 
phase of planetary formation. Our results stress the fact that 
some systems can complete their water transport in a short 
time (~ 6 Myr), in contrast to single star system. Finally, our 
study can give a clear guess on the dynamics and the stability of 
objects moving under the perturbations of a binary star system 
and a gas giant. Indeed, as the depletion of an asteroid belt in 
binaries with small periapsis separations is faster, only a few 
small bodies will remain member of this belt. Thus, it would be 
unlikely to observe such an asteroid belt in such systems. 

It is clear that dynamics in single and binary systems are com¬ 
pletely different as the presence of the secondary together with 
a gas giant would directly have an impact on the dynamics of 
an asteroid belt but also on any plan ets or embryos located in th e 
HZ. Indeed, similarly to the work of lPilat-Lohinger et alJll2008l] . 
secular perturbations and MMRs’ intensity will depend on the 
secondary’s orbital parameters and mass, which will be the pur¬ 
pose of a future study. 
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